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ABSTRACT 

Most of the baryons in the present-day universe are thought to reside in intergalactic space at tem- 
peratures of 10 5-7 K. X-ray emission from these baryons contributes a modest (~10%) fraction of the 
~ 1 keV background whose prominence within the large-scale cosmic web depends on the amount of non- 
gravitational energy injected into intergalactic space by supernovae and AGNs. Here we show that the 
virialized regions of groups and clusters cover over a third of the sky, creating a source-confusion prob- 
lem that may hinder X-ray searches for individual intercluster filaments and contaminate observations 
of distant groups. 

Subject headings: cosmology: diffuse radiation — intergalactic medium — X-rays: general 
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1. INTRODUCTION 

Most of the baryons in the universe remain undetected. 
We believe they exist because primordial nucleosynthesis 
predicts a baryonic matter density amounting to a few per- 
cent of the critical density (p C r), while the baryons associ- 
ated with stars and gas in galaxies make up less than half 
a percent of p cr (e.g., Fukugita, Hogan, & Peebles 1998). 
Within clusters of galaxies the intergalactic baryons are 
obvious because gravitational compression causes them to 
glow prominently in X-ray light. Elsewhere they are much 
more difficult to see. 

Simulations indicate that a large proportion of the uni- 
verse's baryons currently reside outside of clusters in the 
form of diffuse 10 5 ~ 7 K gas associated with groups and 
filaments of galaxies (e.g., Cen & Ostriker 1999; Dave et 
al. 2000). The low surface brightness of this warm gas 
has so far made it very challenging to study. Apart from 
some tantalizing regions of enhanced X-ray surface bright- 
ness (Wang, Connolly, & Brunner 1997; Kull & Bohringer 
1999; Scharf et al. 2000) and a few O VI absorption fea- 
tures (Tripp, Savage, & Jenkins 2000) we have no positive 
detections of it. 

Taken together, the intergalactic baryons lying inside 
and outside of clusters ought to contribute a non-negligible 
fraction of the ~ 1 keV X-ray background. If these in- 
tergalactic gases were heated by purely gravitational pro- 
cesses, they would contribute at least 30% of the observed 
~ 1 keV X-ray background (Pen 1999; Wu, Fabian, & 
Nulsen 1999). However, the point-source contribution at 
these energies is now estimated to be at least 80% (Miyaji, 
Hasinger, & Schmidt 2000; Mushotzky et al. 2000), im- 
plying that less than 20% of this background can be truly 
diffuse. Non-gravitational energy input by supernovae and 
perhaps active galactic nuclei can alleviate this discrep- 
ancy because these energy sources raise the entropy of 
the intergalactic gas, making it harder to compress and 
thereby decreasing its contribution to the X-ray back- 
ground. 



Analyses of the temperature-luminosity relation of clus- 
ters and groups likewise suggest that a significant amount 
of non-gravitational heating has occurred (Evrard & Henry 
1991, Kaiser 1991). The energy injected per particle ap- 
pears to be comparable to the potential depth of a typ- 
ical group, which would severely affect the properties of 
groups (e.g., Ponman, Cannon, & Navarro 1999). Models 
in which early energy injection by supernovae preheats the 
intergalactic medium can plausibly account for this non- 
gravitational heating. The motivation for detecting inter- 
galactic baryons outside of clusters is therefore twofold: 
Detecting this matter would reveal where the baryons im- 
plied by primordial nucleosynthesis have gone, and mea- 
suring the entropy level of these baryons would provide 
a key constraint on supernova energy injection into inter- 
galactic space. 

This paper shows that detecting X-ray emission from in- 
tercluster structures may prove quite challenging because 
the virial radii of groups and clusters out to high redshift 
cover a significant fraction of the sky. Much of the cosmic 
X-ray surface brightness, once point sources are removed, 
will come from this confused patchwork of groups and clus- 
ters, making it difficult to identify true intercluster gas and 
to separate group emission from other virialized sources 
along a common line of sight. 

2. INTERCLUSTER MEAN FREE PATH 

Among the major sources of uncertainty that could com- 
plicate observations of the low surface-brightness struc- 
tures in which intergalactic baryons reside are other col- 
lapsed structures along the same line of sight. The prob- 
ability that any given sight line will encounter a virial- 
ized structure with kT > 0.5 keV turns out to be of or- 
der unity. Here we will estimate that probability using a 
Press-Schechter approach. 

For simplicity, let us assume that the distribution of 
primoridal perturbations is Gaussian with a dispersion 
oc M~ a on mass scale M. The comoving differential 
number density of virialized objects in the universe can 
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then be approximated by dn/dv c = no{v c /vo)~ 1 l a e~ lJ <=l' 2 
(Press & Schechter 1974). In this expression, v c is 
the critical threshold for virialization of a perturba- 
tion at mass scale M in units of standard deviations, 
vq is the value of v c at some fiducial mass scale Mq, 
and n = (2/tt) 1 / 2 (3H^ m /8ttGM ), where H = 
100 h km s _1 Mpc -1 is the current value of the Hubble con- 
stant and Qm is the matter density of the universe in units 
of the critical density. 

An appropriate cross-section for a virialized object is 
, where r2oo is the radius within which the 
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mean overdensity is 200 times the critical density. If we 
define the mass of a virialized object to be the mass within 
r 200 , then cr 200 (M) = n(GM / 100 H 2 ) 2 / 3 , where H is the 
Hubble constant at the redshift of interest. Thus, a mean 
free path A200 between virialized objects more massive 
than M can be defined by 



A 20o( M ) = / O-200 "J"- di/ c 
Jv c (M) av c 

= 2-Pa Q n oV 1 /3a T[l-p^ 2 (M)/2] 



(1) 



where T(a,x) is an incomplete gamma function, p = 
(1 + 3a)/6a, and a = cr 2 oo(M ). 

In order to relate this equation to observations, we 
choose To = 5 keV to be our fiducial temperature 
and apply the low- z mass-temperature relation: kTx = 
(8keV)(M/10 15 h' 1 Mq) 2 ' 3 (Evrard, Metzler, & Navarro 
1996; Horner, Mushotzky, & Scharf 1999). The fiducial 
mass and cross-section then become M = 5x10 14 h~ x Mq 
and 00 = 5.2 h~ 2 Mpc 2 . To find uq, we fit the integral tem- 
perature function j v ( To ^(dn/ dv c ) dv c to the number den- 
sity of > 5 keV clusters, yielding i/q w 2.8 (Donahue & 
Voit 1999). Thus, for a = 1/4, which corresponds to a 
perturbation slope in wavenumber space of = —1.5, we 
arrive at 
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(2) 



In other words, the probability that a line of sight through 
the low-redshift universe will intercept a > 5 keV cluster 
is - 10~ 2 . 

Groups with temperatures of 0.5-2 keV are much more 
common than rich clusters and correspondingly cover a 
much higher percentage of the sky. To estimate that cover- 
ing factor, we extrapolate C/H0X200 down to temperatures 
- 0.5keV assuming awl/4 - 1/6 (n k w -1.5 to -2.0). 
For these cooler objects, we find c/Hq\ 2 oo ~ 0.6 — 1.3. Be- 
cause this quantity is of order unity, the projected virial 
radii of objects with potential depths ;> 0.5 keV ought to 
significantly overlap one another. 

This calculation is, of course, quite simplistic in that 
it fails to account for factors like clustering of virialized 
objects, evolution in their comoving number density, the 
geometry of the universe, and the known tendency for the 
Press-Schechter formula to overestimate the abundance of 
group-size halos. Nevertheless, it illustrates two important 
points: (1) observations of distant groups will be corrupted 
at some level by superpositions with other virialized ob- 
jects, and (2) distinguishing true intercluster baryons from 
those within virialized structures may prove to be quite 
difficult, even under ideal observational conditions. 



3. COVERING FACTOR OF SIMULATED GROUPS AND 
CLUSTERS 

Large-scale cosmological simulations offer a way to as- 
sess the confusion of projected groups more rigorously. 
The largest such simulations performed to date are the 
Hubble Volume simulations of the Virgo Consortium, 
which model comoving cubes 2 — 3 Gpc on a side, suf- 
ficient to reproduce lines of sight stretching most of the 
way through the observable universe (Evrard 1999; Mac- 
Farland et al. 1999; Frenk et al. 2000). Because of the 
magnitude of the computational task, the simulations are 
baryon-free, tracking only the gravitationally-driven be- 
havior of dark-matter particles. Here we will show that 
total solid angle covered by the virial radii of groups and 
clusters resolved in these simulations can approach half of 
the entire sky. 

Hubble Volume simulations have been performed for two 
cosmologies: tCDM (fi M = 1.0, n A = 0.0) and ACDM 
(Am = 0.3, CIa = 0.7). Several "light-cone" catalogs, 
which list the mass and position of each object within the 
past light-cone of a virtual observer located within the 
simulation volume, have been compiled from each simu- 
lation, and the catalogued objects correspond to spher- 
ical regions whose mean mass density is 200 times the 
critical density. Because the masses of individual parti- 
cles in these simulations are w 2 x 10 12 hr 1 Mq, a typi- 
cal group contains < 100 particles. The catalogs list ob- 
jects down to a lower limit of 12 particles, too few to 
model the internal properties of a group but enough to 
qualify as a significant density concentration. An angular 
size for each simulated cluster or group can be computed 
from its catalogued mass and redshift. The coordinate 
distance of such an object is r(z) = J Q 2 c/H(z) dz, where 
H(z) = i?o[^M(l + z) z + SIa] 1 / 2 in a flat universe, and its 
physical diameter is 2r2oo- These quantities combine to 
give an angular diameter of 2r 2 oo(l + z )/ r i z )- 

The most convenient type of Hubble Volume catalog for 
the purpose of studying projection effects is the octant 
catalog which records cluster properties along a light cone 
viewed from a corner of the simulation cube. For ACDM 
and tCDM the limiting redshifts of the octant catalogs 
are z max = 1.45 and 1.25, respectively. Summing the pro- 
jected solid angles of clusters out to the limiting redshift 
and mass of each catalog, we find that the total area cov- 
ered by virialized objects is 46% of the sky in the ACDM 
case and 52% of the sky in the rCDM case. Figure 1 
shows the projected virial radii of clusters and groups more 
massive than 2.7 x 10 13 Mq (12 particles) within a typi- 
cal square degree of the ACDM simulation. Objects with 
< z < 0.5 arc shown in blue, those with 0.5 < z < 1.0 
are in green, and those with 1 < z < 1.45 are in red. Sig- 
nificant overlap is evident in many cases, and accounting 
for that overlap lowers the net sky coverage to 33% for 
ACDM and 35% for tCDM. These figures should be con- 
sidered lower limits for two reasons: Clusters and groups 
more distant than the simulations' redshift limits will pro- 
vide additional coverage, and comparisions of the Hubble 
Volume simulations with smaller-scale ones indicate that 
the Hubble Volume catalogs are 10-20% incomplete at the 
low-mass end (Jenkins et al. 2000). 

Because highly redshifted groups do not contribute 
much flux to the ~ 1 keV background, the apparent tern- 
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peratures of projected virialized objects are also worth in- 
vestigating. Using an evolving mass-temperature relation 
to convert from mass to temperature (Voit 2000) and di- 
viding each temperature by 1 + z, we can assign an appar- 
ent temperature to each cluster and group in the Hubble 
Volume catalogs. Figure 2 shows the same field as Fig- 
ure 1, color-coded by apparent temperature. Objects with 
kT I (1+z) > 1 keV are shown in blue, those with 0.5 keV < 
kT/(l + z) < IkeV in green, 0.5 keV < kT/(l + z) < IkeV 
in red, and kT/(l + z) < 0.25 keV in magenta. The ma- 
jority of the sky coverage in this catalog evidently comes 
from objects with fcT/(l + z) > 0.25 keV. 

4. SURFACE BRIGHTNESS OF VIRIALIZED OBJECTS 

Even though virialized objects cover a large per- 
centage of the sky, much of that solid angle corre- 
sponds to low surface-brightness emission. Current es- 
timates of the mean ~ 1 keV surface brightness at- 
tributable to intergalactic baryons are in the 1 — 4 x 
10 _16 ergcm -2 s _1 arcmin -2 range (Wang & Ye 1996; Wu 
et al. 1999; Pen 1999; Dave et al. 2000). In this sec- 
tion we crudely estimate the surface brightness /200 of a 
virialized object at its virial radius and show that it is 
probably within an order of magnitude of the mean dif- 
fuse surface brightness, implying that these objects blend 
into the background in the neighborhood of their virial 
radii. 

X-ray luminous clusters generally have surface bright- 
ness profiles that are adequately represented by the rela- 
tion I x (b) = I Q [1 + (&/&c) 2 ]~ 3/3+1/2 , where b is the pro- 
jected angular radius from the center of the cluster, b c is 
the angular equivalent of the physical core radius r c , and 
the parameter (3 reflects the behavior of Ix at large radii 
(Cavaliere & Fusco-Femiano 1978). Integrating this ex- 
pression for a typical value of (3 — 2/3 gives the cluster's 
X-ray flux: Fx — 27rfo 2 io- For 6(7200) 3> b c , we then obtain 

J aoo«ra-f^- )(l + ^r 4 , (3) 

87r r 200 \ r 2oaJ 

where Lx is the X-ray luminosity of the cluster. 

To derive a numerical value for /2001 we adopt 
the core radius-luminosity relation of Jones et al. 
(1998), r c = (125/i- 1 kpc)(L x /l-2 x 10 44 h~ 2 ergs" 1 ) ' 2 , 
and assume T2oo T 1 / 2 !!^ 1 . We also construct 
a ROSAT-band luminosity-temperature relation from 
the relations of Markevitch (1998): L x = 2.5 x 
10 43 h~ 2 ergs' 1 (T/2.3keV) c , where < = 2.02 for T > 
2.3 keV and C = 2.64 for T < 2.3 keV. For 
the fiducial temperature of 5 keV, we obtain /200 ~ 
10~ 15 crgcm~ 2 s _1 arcmin -2 (H/H ) 3 {\ + z)~ 4 , very chal- 
lenging for ROSAT but well within the capabilities of 
XMM/Newton and Chandra. Note also that the H 3 
factor owing to the increased high-z density scale tends to 
offset the (1 + z) surface-brightness dimming. 

On group scales, /200 drops by an order of magnitude to 
below the presumptive mean surface brightness of the in- 
tergalactic medium. However, this scaling presumes that 
the luminosity-temperature and core radius-luminosity re- 
lations we have adopted can be extrapolated to <J 1 keV. 
Observations of the Lx-T relation for groups show that 

1 While this paper was being refereed, Croft et al. (2000) released a 
that support our conclusions regarding the difficulty of isolating indn 



this relation in fact steepens at low temperatures and its 
scatter widens (e.g., Helsdon & Ponman 2000). The scal- 
ing of r c is even less secure, given the degeneracy between 
(3 and r c in fits to flattened surface-brightness profiles. 

Combining these surface brightness estimates with the 
Hubble Volume catalogs, we can also assess the distribu- 
tion of X-ray surface brightness owing to virialized ob- 
jects. Figure || illustrates the fractional area fl(> Ix) 
covered by diffuse, virialized gas with a ROSAT-band 
surface-brightness greater than Ix- Splitting the sample 
into objects with masses less than 10 14 h~ x M© and those 
greater than that value shows that the low-mass objects 
dominate the diffuse surface-brightness distribution below 
10 -15 ergs -1 cm~ 2 arcmin~ 2 , raising the possibility that 
the properties of these lower-mass systems might prof- 
itably be studied through statistical analyses of the dif- 
fuse X-ray background, circumventing the difficulties of 
selecting and modelling X-ray groups in an unbiased way. 
In practice, however, such studies will have to contend 
with background counts generated by charged particles, 
which can be comparable to the astrophysical background 
in Chandra and XMM/Newton observations. 

These crude estimates show that while the surface- 
brightness profiles of rich clusters should dominate the 
mean background emission from other diffuse structures 
to at least r2oo> the same is not necessarily true for groups. 
Emission from near the virial radii of groups could be se- 
riously diluted by other emission from unassociated gas 
along the same line of sight, particularly that from other 
groups seen in projection. Numerical simulations account- 
ing for these projection effects will be needed to establish 
the level of this contamination. 1 

5. SUMMARY 

This paper has shown that the projected virial radii of 
groups and clusters cover much of the sky. Estimates 
of this covering factor based on Press-Schechter analy- 
sis place the covering factor of objects p> 0.5 keV near 
unity. Calculations of the group-cluster covering factor 
based on the Hubble Volume simulations corroborate this 
estimate. The virial regions of simulated objects with 
kT/(l + z) > 0.25 cover at least a third of the sky in both 
rCDM and ACDM cosmologies, meaning that projection 
effects could potentially contaminate X-ray observations 
of distant groups. These same projection effects will also 
complicate the search for emission from warm baryons as- 
sociated with intercluster filaments. Untangling the large- 
scale structures seen in deep X-ray images may therefore 
require additional information from optical redshift sur- 
veys that trace the skeleton of dark matter upon which 
intergalactic baryons are draped. 

The authors would like to thank Caleb Scharf and 
Megan Donahue for useful suggestions. AEE is supported 
by AST-9803199 and NASA NAG5-8458. GLB is sup- 
ported by NASA through Hubble Fellowship grant HF- 
01104. 01-98A from the Space Telescope Science Institute, 
which is operated under NASA contract NAS6-26555. 
This research uses data products made available by the 
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rint describing numerical simulations of the diffuse X-ray background 
il X-ray filaments among the many groups that tile the X-ray sky. 
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Fig. 1. — Projected virial radii of simulated clusters and groups in a ACDM cosmology color-coded by redshift: < z < 0.5 in blue, 
0.5 < z < 1 in green, and 1 < z < 1.45 in red. The lower mass limit is 2.7 X 10 13 Mq. The polar (9) and azimuthal {</>) coordinates of the 
objects are in units of degrees. Many virialized objects overlap others in the field. The sum of the solid angles subtended by virialized objects 
is 44% of the full square degree. 
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Fig. 2. — Projected virial radii of simulated clusters and groups in a ACDM cosmology color-coded by apparent temperature: kT/(l + z) > 
1 keV in blue, 0.5 < kT/(l + z) < 1 keV in green, 0.25 < kT/(l + z) < 0.5 keV in red, and 0.25 keV < kT/(l + z) in magenta. The lower mass 
limit is 2.7 X 10 13 Mq, and objects with 0.25 < kT/(l + z) < IkeV dominate the sky coverage. 
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Fig. 3. — Total solid angle f2(> Ix) covered by virialized regions with 0.1-2.4 keV surface brightness greater than Ix- Dotted and dashed 
lines show the contributions to €l(> Ix) from objects with virial masses (M200) above and below 10 14 h Mq, respectively. The sky coverage 
levels off at low flux levels because all surface-brightness profiles are truncated at T200- At intermediate values of Ix, the slope of Q(> Ix) 
corresponds to a /3-model with f) = 2/3, or equivalently, Ix oc t> -3 . 



